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The phase structure of QCD with various number of flavors is studied for Wilson quarks. For the case of 
Nf — 3 we find that the finite temperature deconfining transition is of first order in the chiral limit on an N t = 4 
lattice. Together with our previous results that the deconfining transition in the chiral limit is continuous for 
Nf = 2 and is first order for Nf ~ 6, the order of the transition is found to be consistent with a prediction 
of universality. The case of SU(2) QCD is also studied in the strong coupling limit and the phase structure is 
found to be quite similar to the case of SU(3): There exists a critical number of flavors N F and for Nf > N F the 
confinement is broken even in the strong coupling limit for light quarks. N F = 3 corresponding to 7 for 5(7(3). 
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1. INTRODUCTION 

The deconfining transition of QCD is the last 
phase transition which our Universe has enjoyed 
in the past. To study the nature of this transition 
it is essential to include dynamical quarks. In this 
paper we study the phase structure of finite tem- 
perature lattice QCD with degenerate Np Wilson 
quarks. 

Our model has 4 parameters: Np (the number 
of flavors), (3 (gauge coupling), K (hopping pa- 
rameter) and Nt (lattice size in the temporal di- 
rection). For Nt = 4 - 8 we use a spatial 8 2 x 10 
lattice (the lattice size 10 is doubled for hadrons), 
while for N t = 18 we use 18 2 x 24. The method 
of simulation is similar to that described in . 

1.1. N F — 2 and 6 

In a previous paper |jj we studied the deconfin- 
ing transition Kt for the cases of Np = 2 and 6 on 
an Nt = 4 lattice. Fig. [j] shows our result of the 
phase diagram for Np = 2 near the chiral limit 
Kc- (The result for Np = 3 will be discussed 
later.) Kt is identified by a sudden change in 
the behavior of physical observables such as pla- 
quette and m 2 . Kc is determined by a linear 
extrapolation (in 1/K) of m 2 . in the confining re- 
gion on the lattice with the same N t . We find 
that Nt- (as well as Np-) dependence of Kc is 
small for (3 < 4.5 compared with the difference of 
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Figure 1. Phase diagram for Np = 2 and 3 on an 
N t = 4 lattice. 



Kc determined by m 2 . and m q ; the difference is a 
result of the 0(a) corrections concerning the chi- 
ral behavior of Wilson fermions. The location of 
the crossing point (3ct is identified by monitoring 
iVinv, the number of CG iterations to invert the 
quark matrix, just on the Xc-line 0,01. This can 
be done because of an empirical law |3j that A^ nv 
becomes very large in the confining phase when 
we get close to Kc, while N lnv remains small in 
the deconfining phase even at Kc- 

The nature of the transition at [3ct can be 
studied by measuring observables on the ifc'-line 
We find that, when we decrease (3 toward 
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Figure 2. (a) and 2m q and (b) Mnv history 
for Nf — 3 on Kc-Vme. Jackknife errors for (a) 
are smaller than the symbols. 



Per, Tn^ decreases linearly and is consistent with 
zero at 0ct- This implies that the transition is 
continuous (second order or crossover), in accord 
with the result of an universality argument by 
Rajogopal and Wilczek who predict the second 
order transition for Np =20. 

This result is in clear contrast with the case 
of Nf = 6 where we find that is large near 
and above 0ct — 0.3. We also find a two-state 
signal in the time development of N- lm and other 
observables at = 0.3 fjjj. The first order tran- 
sition observed is again consistent with a predic- 
tion of an universality argument stating that the 
transition is of first order for Np > 3 || . 

2. N F = 3 

Because the effective Nf at the deconfining 
transition temperature in the real world is be- 
tween 2 and 3, it is important to extend the study 
to the case of 3 flavors. Performing simulations 



we find the phase diagram shown in Fig. |l|. As 
expected, Kt shifts to smaller with increasing 
Np from 2 to 3. Fig. || (a) is the result for 
on the Kc line. We find that is large for 
> Act — 3.0. We also find a two-state signal 
at = 3.0 (Fig. ^ (b)). We therefore conclude 
that the deconfining transition is of first order 
for Nf = 3 in the chiral limit, which is in accord 
with the prediction of universality by Pisarski and 
Wilczek [|. 

3. MANY FLAVORS 

3.1. SU{3) case 

Now let us consider the case of much larger Nf ■ 
In our previous paper we have shown that, when 
Np > 7, light quarks are not confining even in 
the strong coupling limit [Q. Extending the study 
to finite we find the following phase structure 
for Np > 70: We have a iVt-independent de- 
confining transition line Kp>, which starts from 
a finite K at = and extends to larger 0, 
In addition, we have a 7V t -dependent deconfining 
transition/crossover line Kt, which reaches the 
i^D-line without crossing the i^-line. 

To understand this phase structure, let us re- 
call the property of the perturbative function. 
Asymptotic freedom is broken for Np > 17 be- 
cause the first coefficient of the function be- 
comes positive. We note that the second coeffi- 
cient changes its sign for Np > 9. This suggests 
the existence of an infrared fixed point (IRFP) for 
16 > N F > N F (N F = 9 up to 2-loop perturba- 
tion theory) . With an IRFP correlation functions 
scale with non-canonical powers at large distances 
and therefore there are no hadronic particles Q|. 
Our MC result jl],| implies that N F = 7. 

3.2. SU(2) case 

As an extension of the SU(3) case, we study an 
SU(2) QCD in the strong coupling limit. We note 
that the SU(2) /3-function has the same charac- 
teristics as 5/7(3) except for that typical num- 
bers for Np are generally smaller than those for 
SU(3): With SU{2), asymptotic freedom is lost 
when N F > 11 and 2-loop result for Np is 6. This 
is natural because the confining force is weaker 
for SU(2). We may therefore expect that N F is 
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Figure 3. Observables in 577(2) QCD for if = 
0.25 at (3 = as a function of Af: (a) plaquette 
and Polyakov loop and (b) and 2m 9 . 



smaller than the case of SU(3). 

Performing a simulation for Np — 8 and N t = 4 
we find that the deconfining transition appears 
at K = 0.2 - 0.21 and that the chiral limit at 
K = 0.25 is in the deconfining phase. Decreas- 
ing Np gradually with fixing K — Re, we find 
that the chiral limit remains in the deconfining 
phase down to Np — 3. When we further de- 
crease Np to 2, iVj nv shows a rapid increase with 
the molecular-dynamics time and finally exceeds 
10 000 in clear contrast with small numbers of 
O(10 2 ) for N F > 3. The chiral limit is therefore in 
the confining phase for Np = 20. Fig. || shows 
our results for physical observables for K = Kc- 
To study the Nt dependence of the transition, we 
simulate on an Nt = 8 lattice for the critical case 
of Np — 3. The stability of the result confirms 
that the deconfining transition we observe is a 
bulk transition. According to our experience with 
SU(3) QCD, we expect that these properties in 



the strong coupling limit hold also for small f3. 
We therefore conclude that N* F = 3 for SU{2) 
QCD. 

4. CONCLUSION 

We studied the phase structure of QCD with 
dynamical Wilson quarks. Extending our previ- 
ous study, we found that the finite temperature 
deconfining transition for Np = 3 is of first order. 
We also studied an SU(2) QCD and found that 
A^-dependence of SU (2) QCD is quite similar to 
that of SU(3), but now N F , above which confine- 
ment is lost for light quarks, is 3 instead of 7 for 
577(3). 

We have studied the phase structure of QCD 
mainly on an Nt = 4 lattice. To study the impli- 
cation to the continuum limit, the system should 
be investigated also at larger (3 on a larger lattice 
|7]|| . We are extending our study in this direc- 
tion. 

The simulations are performed with HITAC 
S820/80 at KEK and with QCDPAX. We thank 
members of KEK for their support and the other 
members of QCDPAX collaboration for their 
help. This project is in part supported by the 
Grant-in-Aid of Ministry of Education, Science 
and Culture (No.62060001 and No.02402003). 
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